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ABSTRACT 

When and how galaxy morphology such as disk and bulge seen in the present-day universe 
emerged is still not clear. In the universe at z > 2, galaxies with various morphology are seen, and 
star-forming galaxies at z ~ 2 show an intrinsic shape of bar-like structure. Then, when did round 
disk structure form? Here we take a simple and straightforward approach to see the epoch when 
a round disk galaxy population emerged by constraining the intrinsic shape statistically based 
on apparent axial ratio distribution of galaxies. We derived the distributions of the apparent 
axial ratios in the rest-frame optical light (~ 5000 A) of star-forming main sequence galaxies at 
2.5 > z > 1.4, 1.4 > z > 0.85, and 0.85 > 2 > 0.5, and found that the apparent axial ratios of 
them show peaky distributions at z > 0.85, while a rather flat distribution at the lower redshift. 
By using a tri-axial model {A > B > C) for the intrinsic shape, we found the best-fit models give 
the peaks of the B/A distribution of 0.81 ± 0.04, 0.84 ± 0.04, and 0.92 ± 0.05 at 2.5 > z > 1.4, 
1A> z > 0.85, and 0.85 > z > 0.5, respectively. The last value is close to the local value of 0.95. 
Thickness {C/A) is ~ 0.25 at all the redshifts and is close to the local value (0.21). The results 
indicate the shape of the star-forming galaxies in the main sequence changes gradually, and the 
round disk is established at around 2 ~ 0.9. Establishment of the round disk may be due to a 
cease of violent interaction of galaxies or a growth of a bulge and/or a super-massive black hole 
resides at the center of a galaxy which dissolves the bar structure. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: structure 


1. Introduction 

When and how did disk galaxies form? At 
z < 1, disk galaxies that can be cla ssifie d with 
the Hubble’s tuning fork exist (e.g., iSchade et al 


metry. Meanwhile, in the universe at z > 3, galax¬ 
ies with various morphology are seen, and many 
of them show the presence of clumps and/or ir¬ 
regularity, though the m orphology is tra c ed in 


2007: Sargent et al.ll20oik they are identified with 

their apparent morphology, and/or their surface 
brightness distribution, and/or empirical indices 
describing such as central condensation and asym¬ 


the rest-frame UV (e.g., iGiavalisco et al. lEael 


—^—II-——a---1|- >— n 1 7 | || - jiiie le sb-iifcime u v ivjicivaT-Lsuu et cll.i 

199o, Abraham et al.|l996| ; |Lilly et al.|l998| ;| Scarlata et al gjteidel et alii 1 996l) . Hence, the emergence of mor 


phology classified with the Hubble’s tuning fork is 
considered to be around at z ~ 1 — 3, and forma¬ 
tion of disk structure would also be expected at 
the epoch. 
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The Wide Field Camera 3 (WFC3) on the Hub¬ 
ble Space Telescope ( HST ) enables us to study 
morphology of galaxies in the rest-frame optica, ! 
light well up to at 2 ~ 3. IConselice et al.l ( 201lh 
studied morphology evolution based on eye inspec¬ 
tion. They found the disk population is increasing 
with redshift gradually and the population is very 
rare at z > 2. By using apparent morphology, s im- 
ilar trend is reported bv ICameron et aD ( 2011 1 for 


galaxie s with a ste l lar m as s range of IQ 10 12 M P 


and bv iTalia et al.l ( 20141) . iMortlock et al. ( 2013 ) 
also studied the morphology evolution with eye in¬ 
spection and a clearer emergence of the disk pop¬ 
ulation at 2 ~ 2 is seen after correcting for the 
statistical morphology misclassification. 

However, from the optical surface brightness 
distribution, th e fraction of the disk galaxies is 
higher at 2 > 2 (Mortlock et al. 20131 1 Bruce et al.l 


120121 Buitrago et al. 20131) . If the Sersic index 


(n) less than 2.5 is regarded to be an indicator of 
disk galaxies, the fraction of the disk population 
is larger at 2 > 2 and is decreasing with redshift, 
though the details of the evolutional trend is not 
so simple; e.g., if irregular galaxies are included 
in the sample, t he disk fraction is m ore constant 
against redshift ( Buitrago et al.ll2013t ). or the frac¬ 
tion is rather constant for less ma s sive galaxies 
(<1O 1O - 5 M 0 ) ( Mortlock et, al. 2013 ). Bruce et al 


( 20121 ) tried to model the two-dimensional surface 
brightness distribution with two components of 
Sersic index of n = 1 and n = 4, and found a 
fraction of the disk (n = 1 component) dominated 
galaxies decreases at 2 ~ 2 rather sharply for stel¬ 
lar mass larger than 10 10 ' 5 Mq. 

Thus the opposite views are obtained in terms 
of the disk evolution. Other parameters charac¬ 
terizing the morphological features such as central 
condensation, asymmetry, Gini, M20, etc seem 
not to be suitable to trace the mor pholo gy of 


Mortlock et al. 2013 : Talia et all 2014 ). 


galaxies in distant universe (IWang et ah 2012; 


On the other hand, studies of internal kine¬ 
matics in star-forming galaxies such as BM/BX 
and star-forming BzK (sBzK) galaxies at 2 ~ 2 
revealed about one third of thcmjshow a cl ear ro - 
tation kinematics ( Forster Sclireiber et al.l[20091) . 
This suggests we witness disk galaxies and they 
may be progenitor of present-day disk galax¬ 
ies, though a fraction of them may evolve into 
the present-day elliptical galaxies through ma¬ 


jor merges, particularly those with large stellar 
mass. Another piece of the suggestion for the 
disk progenitor scenario stems from the cluster¬ 
ing property of sBzK galaxies; sBzK galaxies with 
K < 23 mag are considered to be a progenitor 
of present-day disk galax ies with regard to their 
weak clustering strength (jHavashi et al.ll2007l ). It 
would be reasonable to suppose most of the sBzK 
galaxies in this magnitude range will evolve into 
the present-day disk galaxies, though the brighter 
sBzK galaxies may be a progenitor of elliptical 
galaxies. 

Motivated by the finding of the presence of 
many rotating disk among sBzK galaxies and clus¬ 
tering nature of them, we took a simple and 
straightforward approach to see whether the in¬ 
trinsic shape of faint sBzK galaxies at 2 ~ 2 is 
disk-like or not statistically based on apparent 
axial ratio distribution. We studied the surface 
brightness distribution and the apparent axial ra¬ 
tio in both rest-frame UV and optical lightfor 


sBzK galax i es with K r < 24 mag (Yumc^etjil 


201 lL 120121) . lYuma et al.l (l201ll) and 


Yuma et al 


(12012 ) revealed the surface brightness distribution 


of them is characterized by Sersic index of n ~ 1, 
and the half-light radius and surface stellar mass 
density is similar to those in the present-day disk 
galaxies (Barden et al. 120051 ). However, it was 
turned out that the intrinsic shape of the sBzK 
galaxies is not a round disk (round means intrin¬ 
sic axial ratio of ~ 0.9 — 1.0); their intrinsic shape 
is bar-like or oval (intrinsic axial ratio ~ 0.6 — 0.8) 
with the thickness similar to that of the present- 
day disk galaxies. The intrinsic bar-like structure 
of star-forming gala xies at 2 ~ 2 is also pointed 
out bv lLaw et al.l ( 2012 ). It should be noted here 
that this bar-like structure does not imply a direct 
progenitor of the local barred galaxies as discuss 
later. 

Then, when did round disk galaxies appear af¬ 
ter 2 ~ 2? Emergence of the round disk popula¬ 
tion should reflect the evolution process of the disk 
structure. Hence, revealing evolution of the intrin¬ 
sic shape of the galaxies is expected to give us an 
insight to physical process of the galaxy evolution. 
In this paper, we study the intrinsic shape of main 
sequence galaxies with K§ < 24 mag in the red¬ 
shift range from 2.5 to 0.5 using the same method 
as Yuma et al. (2011) and Yuma et al. (2012) to 
see the redshift that the round disk galaxy pop- 
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illation emerged. In the subsequent section, data 
sources and sample galaxies are summarized. In 
§3, we derive the surface brightness distributions 
and axial ratios of the galaxies. Then, in §4, we 
examine the distributions of the apparent axial ra¬ 
tios and constrain the intrinsic shape of the galax¬ 
ies. Conclusion and discussion are given in §5. We 
adopt a cosmology with H 0 = 70 km s -1 Mpc -1 , 
fl m = 0.3, and Da = 0.7. All magnitudes in this 
paper are given in the AB magnitude system. 


2. Data sources and sample galaxies 


Sample galaxies in this study are collected from 
the Great Observato ries Origins Deep Su rvey- 
South (GOODS-S)Q ( Dickinson et al.l 120031) and 
the Subaru XMM-Ne wton Deep Survey (SXDS) 
( Furusawa et al.ll2008l) . which almost overlaps Ul¬ 
tra Deep Survey (UDS) in UK Infrared Telescop e 
Infrared Deep Sky Survey (UKIDSS) (Lawrence et al 
2007l h The data sou rces a re the same as those 
used by lYuma et al.l ( 20121) and the details are 
described in the paper. Here we briefly summa¬ 
rize them. In GOODS-S region, near infrared 
(NIR) images (J-, H -, and A's-band images) 
were taken from GOODS/ISAAC (Infrared Spec¬ 
trometer And Array Camer a) data release final 
version ( Retzlaff et al.l [201(1) . Five sigma limit¬ 
ing ATs-band magnitude is 24.4 mag. NIR im¬ 
ages taken with WFC3 on HST with F125W 
and F160W filters were obtained from CAN- 
DELS (Cosmic Assembly Near-infra red De e p Ex- 


DtPS /Oosmic Assembly lNe a.r-mtrared De e p Hix- 
tragalactic Legacy Su rvey; iGrogin et al. boilt 


Koekemoer et al data release v0.5. Opti¬ 

cal images were taken from ver. 2.0 data products 
of GOODS HST /ACS(Advanced Camera for Sur- 
veys) treasure program (-B 435 , U 6 60 , * 775 , and z 8 5 o) 


( Giavalisco et al.l 2004) and from ESQ/G OODS 


program (U and R) (iNonino et al. 20091) . Mid- 
infrared (MIR) images at 3.6 /im and 4.5 /.im 
taken with Infrared Array Camera (IRAC) on the 
Spitzer Space Telescope (SST) were from data 
release 1 and 2 obtained in Spitzer Legacy Sci¬ 
ence program. In SXDS/UDS region, NIR im¬ 
ages (J, H , and K$) were taken from data re¬ 
lease ver. 8 of UKIDSS UDS. Five sigma limiting 
magnitude is A's = 24.6 mag. NIR images with 
I7ST/WFC3 with F125W and F160W filters were 


1 http://www.sci.edu/science/goods 

2 http://candels.ucolick.org/index.html 


obtained from CANDELS data release vl.O. Op¬ 
tical images were ta ken from SXDS proje ct ( B , 
V, Rc, i', and z') I Furusawa et al. 2008 ). We 
used HST /ACS images taken with F814W filter 
(* 814 )- MIR 3.6/im and 4.5/mr images were taken 
from Spitzer Public Legacy Survey of UKIDSS 
UDS (SpUDS; PI: J. Dunlop). We also use the 
far-UV and near-UV data tak en from Galaxy 
Evolution Explorer (GALEX: iMartin et al. 2005 ) 
archived image (GR6) and the u-band data of the 
Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS) wide held as supplemental informa¬ 
tion, though the limiting magnitudes are shallow. 

Object detection was made based on AVband 
image by using SExtractor (Beilin & Ar nouts 


19961 ). We made photometry using these data 


for each source detected in K$ images at the 
same position and derived total magnitudes us¬ 
ing SExtractor. For the images taken with ACS 
and IRAC, we applied a correction for the aper¬ 
ture photometry to obtain a total magnitude fol¬ 
lowing the manner by Yuma et al. (2012). Spec¬ 
tral energy distribution (SED) was constructed for 
each object from the photometry in all the bands. 
Then photom etric redshifts were d erived by us¬ 
ing hyperz ( Bolzonella et al.l 120001 ). Compar¬ 
isons with spectroscopic redshifts available showed 
Az/(l + z) of ~ 0.05 and ~ 0.03 in GOODS-S and 
SXDS, respectively. With the photometric red¬ 
shifts, stellar mass, star-formation rate (SFR), and 
color excess ( E(B — U)) were derived for each ob¬ 
j ect th rough SED fitting by using SEDfit program 
( Sawickil[2012 ). The SEDfit em ploys p op ulati on 
synthesis code of BC03( Bruzual fc Charlodl2003 ). 
Salpeter initial mass function with a mass range 
of 0.1 — 100 M 0 and the solar metallicity were as¬ 
sumed. Extinction law used was that by Calzetti 
et al. (2000). More details are described by Yuma 
et al. (2011, 2012) and Yabe et al. (2012, 2014). 

Using these data, we selected galaxies brighter 
than Ks = 24.0 mag. We divided the galaxies in 
three epochs of 2.5 > z > 1.4, 1A> z> 0.85, and 
0.85 > 2 > 0.5 based on the photometric redshifts; 
these epochs are chosen to take a comparable du¬ 
ration (~ 2 Gyr) of the epochs in the cosmic age 
and to keep a reasonable sample size. 

Figure |T] shows distribution of the SFR against 
stella r mass in each epoch. As like in other w orks 
(e.g.. iNoeske et all 12007 : Daddi et, al.ll2007t ). the 
SFRs are larger in more massive galaxies and they 


3 

































































make a sequence referred as main sequence. The 
main sequence shows a gradual cosmological evo¬ 
lution down to 0 ~ 0, where we see disk galaxies 
in the sequence. Most of our sample galaxies at 
2.5 > z > 1.4 are also sBzK galaxies. Galaxies 
in the most massive part of the main sequence 
at the higher redshifts may soon stop star for¬ 
mation and evolve into passive elliptical galax¬ 
ies a ccording to a simple theoretical expectation 
fe.g.. iBouche et ah 2010h . But most of the sam¬ 
ple galaxies occupy less massive part of the main 
sequence and they are expected to reside in the 
main sequence at the present epoch. A branch 
seen above the main sequence shows very high 
SFR at a fixed stellar mass, and they are con¬ 
sidered to be violently formin g stars o ften under 
strong galaxy interaction (e.g.- lDaddi et alJl2007f) . 
Thus we selected star-forming galaxies in the main 
sequence bracketed by thresholds shown in Fig¬ 
ure [1] The stellar mass limit at 5 x 10 9 Mq is 
to avoid the inclusion of Magellanic-type galaxies 
that would tend to show bar-like structure even 
in the local universe. A star-forming galaxy with 
a stellar mass of ~ 5 x 10 9 Mq at z = 2.5 is a 
progenitor of the Milky Way class galaxy (stel¬ 
lar mass of 5 x 10 1 o Mq), if we adopt the abun¬ 
dance m atching in terms of the number density 
( van Dokkum et al. 201 3(1. Ranges of the specific 
SFR (sSFR) we use are 3 x lO” 11 - 2 x 10" 8 
yr" 1 (z = 1.4 - 2.5), 3 x 10" 11 - 1 x 10“ 8 yr" 1 
(z = 0.85 - 1.4), and 2 x 10" 11 - 3 x 10~ 9 yr" 1 
(z = 0.5 — 0.85). SFR limit is just an additional 
threshold; SFR > 2.0M 0 yr -1 ' (z = 1.4 — 2.5), 
> 0.6M 0 yr -1 (z = 0.85 — 1.4), and > 0.2M 0 yr -1 
(z = 0.5 — 0.85). A fraction (typically 10 — 15%) 
of them show multiple-structure in the rest-franre 
optical HST images; if an object detected by 
SExtractor locates within a radius of median 
size (FWHM) of the galaxies in each epoch, we 
regard them as multiple. Such galaxies are not 
adequate to examine the shape and we discarded 
them in this study. 

3. Two-dimensional surface brightness fit¬ 
ting 

We made two-dimensional surface brightness 
distribution fitting to the galaxies selected above. 
The images were taken with HST in the rest- 
frame optical wavelength (~ 0.5^m): F160W(iii6o) 
for 2.5 > 2 > 1.4, F125W( J 125 ) for 1.4 > z > 0.85, 


and F85OLP(0 85O ) in GOODS-S and F814W(i S i 4 ) 
in SXDS for 0.85 > z > 0.5. Magnitude ranges for 
the sample galaxies in three epochs were virtually 
H i 6 o < 24 mag and Ji 25 < 24 mag, and 2 850 < 25 
mag in GOODS-S and isi 4 < 24 mag in SXDS. 
Pixel scales after drizzle were 0."06, 0."06, and 
0."03 in # 160 , di 25 , and z 85 o/i 8 i 4 , respectively. 

The radial distribution is fitted with Sersic pro¬ 
file: I(r) = I e exp[—ftn{(r/r e ) 1 / n — 1}], where r e 
is the half light radius, I e is the surface bright¬ 
ness at r e , and n characterizes the shape of the 
profile (Sersic index). Two dimensional surface 
brightn ess fitting was ma de by using GALFIT ver¬ 
sion 3 ( Peng et al ] m3). Initial guess parame¬ 
ters on the central position, axial ratio, position 
angle, and r e are taken from the output parame¬ 
ters of SExtractor. Initial Sersic index was set 
to be 1.5, which does not affect the final results 
significantly. Point spread function in each image 
was constructed from 5 ~ 20 isolated unsaturated 
stars. 

Figure [2] shows example of the profile fitting. 
As seen in the figure, the surface brightness pro¬ 
files extend to ~ 1 arcsec in radius correspond¬ 
ing to ~ 8 kpc at z = 1 ~ 2 and ~ 6 kpc at 
z = 0.5. The surface brightness we use to fit the 
model reaches 23.5 ~ 25 mag arcsec^ 2 in the rest- 
frame V band after correcting for the cosmolog¬ 
ical dimming. This level of the surface bright¬ 
ness can be regarded as the outskirt of galaxies, 
and is deep enough to confront with the shape 
of local disk galaxies. Considering the fading of 
galaxy luminosity to the present epoch, the value 
woul d be ~ 1 mag arcsec -2 fainter at z ~ 0 
(e.g., Brinchmann et al. 1998; Miller et al ] mil). 
We obtained r e , n, and apparent axial ratio (b/a, 
where a is a major axis and b is a minor axis) for 
each galaxy. About 60-70% of the sample galax¬ 
ies show n = 0.5 — 2.5. We then examined the 
distribution of r e and stellar mass of the sample 
galaxies with n = 0.5 — 2.5. In each epoch, the 
distribution of r e against stella r mass is similar to 
that for local disk galaxies ( Barden et al.l 2005 1. 

Accuracy for the Sersic index as well as for 
the axial ratio was examined by putting artifi¬ 
cial objects into the images used for the analy¬ 
sis. The artificial objects were created by using 
ARTDATA/MKOBJECTS in IRAF. Mock galax¬ 
ies with Sersic index of 1 or 4 were generated in 
both GOODS-S and SXDS with half-light radius 
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0.5<z<0.85 


0.85<z<1.4 


1.4<z<2.5 



Stellar Mass (Msun) Stellar Mass (Msun) 


Stellar Mass (Msun) 


Fig. 1.— Star formation rate versus stellar mass of the sample galaxies at three epochs: From right to left, 
2.5 > 2 > 1.4, 1.4 > z > 0.85, and 0.85 > 2 > 0.5. Large red points are sample in this study and solid lines 
in each panel show the selection criteria 

(see text). 


of 0."1 — l."0 in a magnitude range of 20 — 25 
mag. The ranges of input parameters are cho¬ 
sen to mimic those of the sample galaxies. Ax¬ 
ial ratio and position angle are randomly chosen 
and the galaxies are put into H i 6 o, J 125 , zssth and 
*814 images randomly. As we did for the real ob¬ 
jects, we first used SExtractor to derive posi¬ 
tion, position angle, axial ratio, and FWHM. We 
then performed GALFIT with these parameters 
as well as n = 1.5 as the initial guess. Results 
for iLi6o(GOODS-S, SXDS/UDS) were presented 
by Yuma et al. (2012). Accuracy for the Sersic 
index in ffwo inrage is good enough (< 10 %) 
in this magnitude range in case of n = 1. Even 
if the intrinsic n is 4, the misclassification into 
n = 1 disk is negligible, thou gh the accuracy is 


slight ly worse in z 85 o image (|Yuma et al.1 12011 


20121) . We show the accuracy of the axial ratio de¬ 


rived for J 125 (SXDS/UDS), zs 5 o(GOODS-S), and 
* 8 i 4 (SXDS/UDS) in case of n = 1 in Figured] In 
all the panels, the obtained axial ratios agree very 
well with those of input value in the bright magni¬ 
tude range, and the accuracy is getting worse for 
the fainter objects. However, at the faintest mag¬ 
nitude level of the sample galaxies in each band, 
the accuracy of the obtained axial ratio is less 
than 5% on average, which is smaller than the bin 
size of the histograms of the axial ratio distribu¬ 
tion described below. The accuracy of axial ratio 
in Hifin(G OODS-S, SXDS/UDS) is presented by 


Yuma et al. ( 2012 ). and it is comparable to those 


in other bands. Detection rate in ivs-band image 
against the apparent axial ratio was also exam¬ 


ined and the rate is almost constant in the range 
of 6/a = 0.2 and 1.0 ( Yuma et al. 2011 ). 


• % 

SXD51_475B 

z photo = 0.6700 Fte=2.08[kpc] (1 = 1.43 b/a=0,58 



5XOS3.117457 

ip ho to = 1.0400 Re=3.02[kpc] h = 1.98 b/a=Q.57 


♦ % 


5XDS1 2027 

Z_phOtO = 2.4300 R«=2.30[kpC] 11 = 1.75 b/0=0.49 



Fig. 2.— Example of the surface brightness dis¬ 
tribution fitting at 0.85 > 2 > 0.5(upper), 1.4 > 
2 : > 0.85(middle), and 2.5 > 2 > 1.4 (lower). The 
images are taken with HST * 314 (upper), WFC3 
J 125 (middle), and WFC3 Hiq 0 (lower). In each 
row, from left to right, original image, best-fit 
model image, residual, and radial surface bright¬ 
ness profiles along the major axis are shown. In 
the rightmost panel, the observed surface bright¬ 
ness and the best-fitted profile smoothed with 
point spread function is shown with circles with 
error and solid curve, respectively. Point spread 
function is shown with thin solid curve. 
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4. Axial ratio distribution 


4.1. Distribution of apparent axial ratios 


In order to see whether they have round disk 
structure, we take a simple and straightforward 
approach; i.e., we constrain the intrinsic shape of 
these galaxies statistically based on distribution 
of apparent axial ratios of the sample galaxies at 
three epochs. Figure [I] (upper row) shows the 
apparent axial ratio distributions of the sample 
galaxies with n = 0.5 —2.5 among galaxies selected 
in §2 If the intrinsic shape of disk galaxies is an 
infinitesimally thin circular disk, the apparent ax¬ 
ial ratio distribution is exactly flat. In fact, the 
distribution for ~ 300,000 local disk galaxies is 
al most flat except f o r aro und at b/a ~ 0 and ~ 1 
( Padilla and Straussl 2008b . since real disks have 
a finite thickness and are not perfectly round. As 
seen in the right panel (2.5 > z > 1.4) of the upper 
row of Figure HI the axial ratio distribution peaks 
at ~ 0.45 and sharply decreases to both sides. In 
the middle panel (1.4 > z > 0.85), the distribution 
is slightly flatter, but still shows a slight peaky dis¬ 
tribution. The histograms suggest these galaxies 
do not have a round disk structure. Meanwhile in 
the left panel (0.85 > z > 0.5), the distribution is 
rather flat and is close to that for the local disk 
galaxies, suggesting they have the intrinsic shape 
of round disk. 


SXDS ^ . 

1.0 

GOODS-S Djsk # 

1.0 

SXDS Dilk . 


j? 0-5 


0.5 


-•-*+i+j 

f 0.0 

. 

f 0.0 



£ -0.5 


£ -0.5 



-1.0 


-1.0 



21 22 23 24 25 20 21 22 23 24 25 26 27 21 22 23 24 25 

Input WFC3/F125W Mag Input ACS/F850LP Mag Input ACS/F814W Mag 


Fig. 3.— Accuracy of the obtained axial ra¬ 
tio against magnitude of input artificial object 
in Ji 25 (left), £ 850 (middle), and * 814 (right). Error 
bars represent Ict dispersion of the distribution. 


4.2. Intrinsic axial ratio distribution 

To constrain the intrinsic shape quantitatively, 
we employ a tri-axial model with axial lengths 
of A > B > C following the method by Ryden 
(2004). The face-on ellipticity (e = 1 — B/A) is as¬ 
sumed to be described by log-normal distribution 
with mean of [j,(= In e) and dispersion of ct, while 


the edge-on thickness (C/A) is assumed to be de¬ 
scribed with Gaussian distribution with mean of 
/x 7 and dispersion of ct 7 . With this parameter set, 
the expected apparent axial ratio distribution seen 
from random viewing angles can be calculated and 
can be compared with the observed apparent axial 
ratio distribution. 


The model axial ratio distributions were cal¬ 
culated in a range of /j = —3.95 ~ —0.05, ct = 
0.2 ~ 2.0, /z 7 = 0.1 - 0.98, and ct 7 = 0.01 - 0.35, 
with a step of A /i = 0.15, Act = 0.15, Ap 7 = 0.02, 
and Act 7 = 0.02. With these parameter sets, we 
searched for the best-fit parameter set using y 2 
method with Ict statistics by Gehrels (1986). To 
estimate the uncertainty of the best-fit parame¬ 
ters, we made Monte Carlo realization of the ob¬ 
served histogram of the axial ratio 1000 times and 
derived the best-fit parameter set for each realiza¬ 
tion. We took 68 % confidence level as the uncer¬ 
tainty of the model parameter. The peak value 
of the distribution of intrinsic B/A is given by 
1—exp(/i — ct 2 ). The uncertainty on the peak B/A 
was calculated based on the propagation of errors 
with the uncertainty of fi and ct mentioned above. 


The best-fit models are shown with solid curves 
in Figure [4] (upper row), and the best-fitted pa¬ 
rameters as well as the uncertainties are shown 
in Table 1. Resulting distributions of intrinsic 
axial ratios at the epochs are shown in Figure 
U (bottom row). The best-fit model B/A peaks 
at 0.81 ± 0.04,0.84 ± 0.04, and 0.92 ± 0.05 at 
2.5 > £ > 1.4,1.4 > £ > 0.85, and 0.85 > £ > 0.5, 
respectively, and the last value of B/A is very 
c lose to that for the local disk galaxies (0.94) 
( Padilla and Straussl 12008 ). The mean C/A is 
~ 0.25 at all the epochs an d is close to that of 
the local disk galaxies (0.21) ( Padilla and Straussl 
2008i h The intrinsic shape of the star-forming 


main sequence galaxies at 2.5 > z > 0.85 is bar¬ 
like or oval, while they have round disk structure 
at £ < 0.85, i.e., the emergence of the round disk 
population is £ ~ 0.9, though the precise redshift 
determination is difficult due to the the gradual 
cosmological evolution and small sample size. 


4.3. Robustness of the results 

In order to examine the robustness of the re¬ 
sults, we constructed samples by changing selec¬ 
tion criteria slightly with respect to Ag-band mag¬ 
nitude, stellar mass, sSFR, * 814 /^ 850 /^ 125/^160 
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Axial ratio (b/a) Axial ratio (b/a) Axial ratio (b/a) 


Becoming round at z<0.85 



0.4 0.5 0.6 0.7 0.8 0.9 1.00.0 0.2 0.4 0.6 0.8 1.0 

Intrinsic face-on axial ratio (B/A) Intrinsic edge-on axial ratio (C/A) 


Fig. 4.— Normalized histograms of the observed axial ratios at each epoch (upper panels) :From right to left, 
2.5 > z > 1.4, 1.4 > z > 0.85, and 0.85 > z > 0.5. Numbers in vertical axis are multiplied by 10. Number 
in the parenthesis at the top of each panel refers to the sample size. Best-fit model distributions are shown 
with solid curves. Axial distribution of local disk galaxies from the Sloan Digital Sky Survey Data Release 
6 by Padilla & Strauss (2008) is also shown with dotted curve (left panel). Distribution of best-fit intrinsic 
axial ratios of B/A and C/A (bottom panels). Red solid, green dot dashed, and blue dashed curves refer 
to the intrinsic model axial distribution at 2.5 > z > 1.4, 1.4 > 2 > 0.85, and 0.8 5 > z > 0.5, respectively. 
Black dotted curves show the distribution for local disk galaxies ( Padilla and Strauss! 12008 ). Y-axis is an 
arbitrary unit. Peaks for the B/A distribution are shown by tick marks with uncertainty above the curves. 


Table 1: Best-fit parameter sets for the tri-axial model. 
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z = 0.5 - 0.85 
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0.50lg;lg 

0 26 +0 02 
u ' zu - 0.16 

0 QQ+ 0 - 01 

u.uy_ 0 .08 


‘Padilla & Strauss (2008) 


magnitude, and Sersic index, and made the same 
fitting to these samples. 

If we take samples in each epoch with 0.5 mag 
brighter or fainter in A'g band, virtually all the 
best-fit parameters are within lcr error. The peak 
values of the intrinsic axial distribution {B/A) 


agree with the best-fit value within 1 cr error. If 
we take a massive sample (stellar mass larger than 
1.0 x 10 10 Mq), the results mostly agree within 1 
a error. If we include less massive galaxies (stellar 
mass larger than 3 x 10 9 Mq), peak values of B/A 
agree within the error, but the resulting parame- 
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ters for the highest-redshift bin change more than 
1 a. It should be noted that even in local uni¬ 
verse such less massive galaxies often show bar- 
like structure as like Magellanic irregular galax¬ 
ies. We also test the dependence on sSFR by tak¬ 
ing 3 x 10 -11 yr -1 < sSFR < 2 x 10 -8 yr -1 at 
z = 2.h- 1.4, 3 x KT 11 yr” 1 < sSFR < 5 x 10" 9 
yr" 1 at 2 = 1.4-0.85, and 3x 10“ n yr” 1 < sSFR 
< 7 x 10 -9 yr -1 at z = 0.85 — 0.5. Almost all the 
obtained best-fit parameters as well as the peak 
values of B/A agree within ler error. If we take 
the threshold SFR 2-3 times larger, the resulting 
parameters are almost within 1 a. The axial ra¬ 
tios are derived by using * 814 , £ 850 , J 125 , and _ffi6o 
images. More accurate axial ratios are expected 
to be obtained with the brighter sample. Hence 
we take about 0.5 mag brighter samples (*si 4 < 

23.5 mag, £sso< 24.5 mag, Ji 25 < 23.5 mag, and 
Hieo< 23.5 mag) from the I\$ < 24.0 mag sam¬ 
ple. The resulting peak axial ratios are within or 
marginally within 1 er uncertainty, and the conclu¬ 
sion in this study does not change. We also make 
samples with Sersic index of 0.5 —1.5 and 0.5 — 3.5, 
and find again the conclusion does not change. 


Dust obscuration may affect the axial ratio dis¬ 
tribution; the redder disk galaxies tend to show 
the smaller axial ratios in t he loc al univ erse as we ll 
as at £ = 0.6 — 0. 9 (e.g., Padilla and Strauss! 2008 : 
Patel et al . 2012). This trend is considered to be 


due to the larger dust obscuration in more inclined 
galaxies. In order to see the effect, we divided the 
sample into two subsamples with smaller or larger 
color excess ( E(B — V)) at each epoch. The color 
excess was derived from the SED fitting mentioned 
in §2, and we divided the sample at around the 
median E{B — V) of 0.2 or 0.25 mag. Resulting 
histograms of the axial ratios show that a slight 
excess at smaller axial ratio (6/a ~ 0.2 — 0.4) is 
seen at £ = 0.5 ~ 0.85, but such a clear trend 
is not seen at the higher redshifts. We made the 
model fitting to the subsamples. About half of the 
best fit /r and a values agree with those of the total 
sample within the uncertainty for the total sample. 
However, it should be noted that the uncertainty 
for the subsamples is large due to the smaller sam¬ 
ple size. We also derived the peak B/A with the 
subsamples. The peak B/A values in the smaller 
color excess (larger color excess) are 0.81 (0.93) at 

2.5 > £ > 1.4, 0.91 (0.78) at 1A > z > 0.85, and 
0.99 (0.97) at 0.85 > £ > 0.5. Thus the increas¬ 


ing trend of B/A with decreasing z is still seen, 
except for the B/A for the larger color excess at 

2.5 > £ > 1.4. 


5. Conclusion and discussion 

We constrain the intrinsic shape of the star¬ 
forming main sequence galaxies at £ = 2.5 to 
0.5 based on the distributions of the apparent 
axial ratios by employing the tri-axial model of 
A > B > C. At £ > 0.9, the intrinsic shape of 
them is bar-like or oval with B/A ~ 0.8, while at 
£ < 0.9, B/A is 0.92 that is close to the value for 
the local galaxies (0.95). Hence, the emergence 
of the round disk can be regarded to be £ ~ 0.9, 
though the precise redshift determination is diffi¬ 
cult due to the gradual evolution and the small 
sample size. The thickness {C/A) is about 0.25 
and does not show significant cosmological evo¬ 
lution down to £ ~ 0. Similar result has been 
obtained very recently by adopting a slightly dif¬ 
ferent model for the intrinsic shape; the triaxiality 
of star-forming galaxies with stellar mass larger 
than 1 0 9 - 5 Min is large at z > 1, but small at z < 1 
( van der Wei et ahl 2014) . It is worth noting the 
bar-like structure we see at the high redshifts is 
unlikely to be a direct progenitor of bar struc¬ 
ture in the present-day barred galaxies, because 
the apparent axial ratios are obtained mainly from 
outside beyond the bar region in the present-day 
barred galaxies; the bar region in a present-day 
barred galax y is brighter than ~ 22 m ag arcsec -2 
in V band ( Ohta et al. 1990L 2007 1. Further¬ 
more, large barre d ga l axies are disappear ing as 
redshift increases ( van den Bergh et al. 1996l f: Its 
fraction among disk galaxies is only 10-20% at 
£ = 0.8 ~ 1.0 ( Sheth et al. 2008 : Melvin et al 


2014) . 

At £ > 3, most of star-forming galaxies 
show Sersic index of ~ 1 in the rest-frame UV 


(|R avindr anath e t, al.l 1200611 and in the rest-frame 
optical (Akiv ama et ah 1200811 . However, the 
distribution of the apparent axial ratios is not 
flat, suggest ing the el onga ted intrinsic structure 
(1R avindranath et al.l 120001 1. The results obtained 


here reveal that although the disk-like structure 
is seen among star-forming galaxies at £ > 1, 
its intrinsic shape is bar-like or oval. Physi¬ 
cal mechanism to make such structure is not 
clear, but they may be bar-like due to galaxy 









































interaction and/or m inor merge (jNoguchil 11987 ; 
iBerentzen et al. [2004) . Or the baryonic matter 
might have been already highly concentrated in a 
dark matter halo, resultin g into the bar instability 
( Ostriker fc Peebleslll973lf . 

Then how did they become round? Numer¬ 
ical simulations show a central mass condensa¬ 
tion (either of a bulge and/or a super-massive 
black hole (SMBH)) with a mass fraction of 1- 
10 % of a disk dissolves the bar-like s tructure 
( Athanassoula et abl 12005 ; Hozumil 2012 ). Thus 
the emergence of round disk population may 
be related to the growth of a central bulge 
and/or a SM BH. Acco rding to the simulations 
by [Athanassoula et al.l ( 2005h . the bar-like struc¬ 
ture gradually evolves into a round disk with a 
time scale of a few Gyr, which is comparable to 
the elapsed time from z = 1.5 to 0.9 (~ 2 Gyr). 
However, eye inspection of our sample galaxies 
at z = 0.85 — 0.5 suggests that bulge dominated 
galaxies are not the majority. The two-component 
model fit to discriminate bulge and disk in a galaxy 
is not straightforward and it is beyond our scope in 
this s tudy. AGN number density peaks at z = 1 — 2 
le.g. JUeda et al.1120031 : lHasinger et al.ll2005l) . sug¬ 
gesting the rapid growth of SMBHs during the 
epoch. The growth of a SMBH may also play a 
role for dissolving the bar structure, although the 
mass fraction may not be large enough. Alterna¬ 
tively, the gradual cease of galaxy interaction may 
make a round disk. If this is the case, the epoch 
of the round disk population at * ~ 0.9 would in¬ 
dicate the major galaxy interaction mostly ceased 
around at this redshift and the galaxy morphology 
is rather stable after this. Further studies of link 
between intrinsic shape of disk and the growth of 
a bulge/SMBH as well as galaxy interaction are 
desirable to disclose the evolution of disk galaxies. 
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